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摘要 : 对 织 氨 原子 兰 姆 位 移 的 测量 是 测量 质子 电荷 半径 的 最 有 效 手段 之 一 ， 而 对 于 绿 氨 原子 兰 姆 位 移 的 理论 
计算 则 是 从 兰 姆 位 移 的 实验 结果 中 抽取 出 质子 电荷 半径 的 必要 途径 。 本 文 计算 了 缪 氨 原 子 的 兰 姆 位 移 ， 考 虑 
了 相对 论 效应 对 组 氢 原 子 兰 姆 位 移 的 影响 。 由 真空 极 化 效应 导致 的 兰 姆 位 移 由 Uehling 等 效 势 能 来 表征 , BA 
原子 的 相对 论 效应 通过 Breit-Pauhli 哈密 顿 量 来 表达 ; 本 文通 过 精确 的 少 体 计算 方法 -高 斯 展开 法 ， 计 算 了 织 
ARF WH 2S 与 2Pi 的 能 级 ， 发 现 考虑 了 相对 论 修正 后 的 兰 姆 位 移 比 不 考虑 相对 论 修正 的 兰 姆 位 移 大 0.09 
meV， 且 更 接近 实验 结果 。 
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引言 能 级 差 称 为 兰 姆 位 移 。 后 来 随 着 量子 电动 力学 (QED) 的 


2 发展， 兰 姆 位 移 才 得 到 解释 中， 其 被 解释 为 电子 与 电磁 
FUR J <Æ H 然 界 最 轻 的 原 J; E 个 由 电 J 和 质 J 场 的 真空 波动 相互 作用 而 产生 的 现象 。 
组 成 的 系统 。 尽 管 是 构成 最 简单 的 原子 ， 但 是 对 其 原子 ee Aeon ae 
光谱 的 探索 却 直接 派生 出 了 现代 物理 学 的 两 大 支柱 之 一 人 所 原子 是 一 个 由 负 缪 子 和 质子 组 成 的 系统 ， 其 寿 


的 量子 力学 。 在 玻 尔 理论 成 功 描述 了 氢 原 子 系统 GER an EF 8 2 E T, = 2.19703(4) x 106s RE. AN 
能 级 ) ja, DADE SEIT ET DUS RR Bh ite KA EHS a Y 200 倍 ， 所 以 其 基态 的 结 


> : ‘ Š ss BATHE ASct BE AK 位 FR 
原子 波 函 数 。 然而 ， 在 后 续 实验 中 发 现 毛 原子 光谱 中 红 合 能 约 为 氨 原 子 基 态 结合 能 的 200 倍 ， 27 keV 左右 ， 对 应 


色谱 线 是 一 个 紧密 间隔 的 双 峰 ， 这 种 分 裂 被 称 为 精细 结 的 缪 氧 原子 的 玻 尔 半径 也 大 约 为 氧 原子 的 1200， 即 265 


构 是 电子 自 旋 和 轨道 角 动 量 耦 合 导致 的 这 是 电 J = 旋 fm 左右 。 同时 ， 缘 子 波 函 数 受 iF 的 影响 大 约 是 电 | 波 


存在 的 早期 实验 证 据 。 为 了 理解 该 现象 ， 狄 拉克 于 1928 。 图 数 受 质 了 影响 的 Cmy/m。) x 10 EM. 不 同 于 氧 原 了 ， 


年 提出 了 Dirac 方程 。Dirac FRETS «TS SUR EINEM Rh, DR TEXT He 


AE, WATER 1 有 2 粒子 的 运动 ， 同 时 考虑 了 相对 论 。 FRR PRR BET, UC SURE IERS 
效应 和 量子 力学 效应 ， 成 功 解释 了 电子 自 旋 现 象 。 但 是 RTS, RET ASB ALS MA 


=I 


` 、 i= fobs 个 全 74 = 机 了 [7-9] SEI IEA 更 o 得 一 ` 
Dirac 方程 计算 出 的 氢 原 子 25, 态 和 2P 态 是 简 并 的 ， 子 结 构 与 电荷 分 布 “” 等 效应 的 影响 更 大 。 值 得 一 提 的 
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NIGER SEU, OTE 2S) 态 和 2P， 之 间 的 真空 极 化 效应 夭 担 SETA PS See A 99% AY 
贡献 ， 而 与 之 有 很 大 不 同 的 是 ， 真 空 极 化 效应 对 氢 原 子 
兰 姆 位 移 的 贡献 比例 很 小 外。 
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除 此 之 外 ， 显 著 影 响 缪 氨 原 子 兰 姆 位 移 的 是 质子 的 电荷 
半径 。 因 此 ， 结 合理 论 上 给 出 的 缪 子 氨 原子 兰 姆 位 移 与 
质子 电荷 半径 的 关系 上 9， 可 以 通过 实验 测量 缪 氧 原子 
289 —2P p 兰 姆 位 移 来 准确 的 给 出 质子 的 电荷 半径 。 9 
在 相关 的 缪 子 氨 原子 的 实验 中 ， 其 给 出 的 质子 电荷 半径 
为 0.84184 土 0.00067 fm"), 但 是 , 另 一 方面 , 质子 的 电 
荷 半径 还 有 另外 两 种 较为 常见 的 测量 方法 ， 即 利用 电子 
质子 弹性 散射 己 -23] 或 氨 原 子 兰 姆 位 移 必 29。 然而 , 这 两 
种 测量 方法 在 不 同 的 实验 中 给 出 质子 电荷 半径 呈现 出 较 
大 的 差异 ; 2010 的 电子 质子 散射 实验 与 2018 EM ZUR 
子 能 谱 测量 实验 户 29 给 出 的 质子 电荷 半径 都 在 0.87 到 
0.88 fm 之 间 ， 而 2019 年 熊 伟 志 等 人 完成 的 电子 质子 散 
射 实验 P 与 2019 年 Bezginov EA CS 完成 的 氧 原子 能 
谱 实 验 给 出 的 质子 电荷 半径 分 别 为 0.831 + 0.019 fm 与 
0.833 + 0.010 fm, 尽管 最 新 的 电子 质子 实验 给 出 的 质子 
电荷 半径 更 为 接近 缪 子 所 原子 给 出 的 质子 电荷 半径 值 ， 
但 是 关于 质子 电荷 半径 的 一 系列 理论 与 实验 研究 仍 需 进 
行 。 

考虑 到 缘 毛 原子 兰 姆 位 移 在 质子 电荷 半径 测量 中 的 
重要 性 ， 也 为 了 更 好 的 研究 缘 子 氢 原 子 的 能 谱 ， 精 确 计 
BOT ART N= We AER, RMB IT 
HPE ESETRE SURF 2S SA 2P 态 
AER, WEHE, HACA HES TZAT 
质子 约 化 质量 的 影响 ， 真 空 极 化 效应 则 通过 Uehling 等 
效 势 能 来 表征 。 通 过 高 斯 展开 法 R11， 计算 上 述 两 种 效应 
WEF ZUR F 2S. 和 2P 能 级 的 影响 ， 并 将 结果 与 一 
阶 微 扰 的 计算 结果 进行 对 比 。 

本 文具 体内 容 如 下 ， 第 二 节 给 出 了 描述 氧 原子 系统 
的 薛 定 订 方 程 ， 非 相对 论 近 似 的 狄 拉克 方程 和 描述 缪 子 
氨 原 子 系统 的 相对 论 性 方程 ， 最 后 给 出 了 缪 子 氧 原子 的 
兰 姆 位 移 的 等 效 势 。 第 三 节 详 细 介 绍 了 高 斯 展开 法 以 及 
具体 计算 中 所 需 的 一 系列 矩阵 元 的 推导 。 第 四 节 是 计算 
结果 ， 首 先 展示 了 和 氢 原 子 的 计算 结果 ， 虽 在 测试 高 斯 展 
开 法 的 精确 度 ; 然后 给 出 了 缪 氨 原 子 兰 姆 位 移 相 关 的 一 
系列 结果 。 最 后 一 节 给 出 了 本 文 的 结论 。 
1 氧 原 子 与 缪 氧 原子 系统 哈密 顿 量 
1.1 氢 原 子 系统 

氧 原子 是 一 个 由 电子 和 质子 组 成 的 二 体系 统 ， 电 子 
HEN m, = 0.5110 MeV, 质子 质量 约 为 电子 质量 的 1800 
倍 ， 远 远大 于 电子 质量 。 所 以 氨 原 子 系统 的 薛 定 谓 方程 


可 以 简化 成 电子 在 质子 的 库伦 引力 势 场 中 的 单 体 运动 方 
程 ， 此 时 握 原 子 的 苹 定 请 方程 如 下 : 


2 
(Z-2)r=er (1) 


Y 即 为 电子 的 波 函 数 , m 为 电子 与 质子 的 约 化 质量 (m x 
m,) fo} 
在 氧 原子 中 ， 非 相对 论 极 限 的 狄 拉克 方程 如 下 : 
| P a pt e? 


2m r 8m 2me2r3 


(s -1) 

(2) 
zhe? / 

+ 下 ao Y=FE Y. 

其 中 除了 和 方程 (1) 中 共有 的 动能 项 和 库伦 势 之 外 ， 另 
外 包括 三 项 相对 论 修正 ， 分 别 为 动能 的 一 阶 相对 论 修正 ， 
自 旋 轨 道 耦合 项 以 及 Darwin 项 。 自 旋 轨 道 耦合 项 只 对 
轨道 角 动 量 1 #0 的 态 有 修正 ， 而 Darwin 项， 又 称 作 接 
触 势 ， 只 对 1 = 0 的 态 有 修正 。 


12 BIRTAR 


BART AB TAM a AM, Al DAR AE BUR 
FRRAPN ES RRNRET. MART E — NEMA 
电荷 ， 旋 为 ;, 这 些 性 质 与 电子 是 相同 的 。 不 同 之 处 在 
于 u 子 质量 为 105.6584 MeV， 约 为 电子 质量 的 200 倍 ， 
XG RE AUR FARR Alm, ca) 要 比 氧 原 子 小 大 
约 200 žo 

NBARTAS, WR AT iA 
子 质量 (m, > m,) 利用 Dirac FEHR AUF ERE 
不 恰当 的 , JABS is SAS hee NR, BART 
ARCEN ERRE RETR FWA IBLE, A 
此 应 该 使 用 两 体 Breit-Pauli GWE), KREN] 
Breit-Pauli 哈密 顿 量 的 部 分 项 ， 具 体 如 下 : 


A 
+ 登 x( SPs 54d ae G) 
My m 


其 中 , m, 是 继 子 和 质子 的 约 化 质量 , m, 是 缪 子 质量 , m, 
是 质子 质量 ，ci CAN Ae. MAE (2) 相 比 ， 其 额 
外 考虑 了 质子 质量 与 质子 半径 带 来 的 修正 。 


1.3 兰 姆 位 移 


兰 姆 位 移 是 2S — 2Pp 的 能 级 差 ， 对 于 氢 原 子 系 
统 ， 对 兰 姆 位 移 的 最 大 贡献 是 电子 的 真空 自 能 ;而 对 于 
BIRTAR, 单 圈 真 空 极 化 约 贡献 了 99% 的 兰 姆 位 移 ， 
本 文 使 用 Yukawa 势 作 为 单 圈 真 空 极 化 的 等 效 势能 中 : 


Vi(r) = Ae Pr. (4) 


其 中 1 = 9.3137 x 10-3 fm™!, A = 2.9533 x 107? MeV. 


‘2. 原子 核 物 理 评 论 第 卷 
= # bes 2: 
2 高 斯 展开 法 Spn = (i aS 51m) 
为 J TR APR DN AS ERE ZE TS 7 FS: e? gq _ 1)! 3 (14) 
R2 mma an e) Nar 
ev +V(r)- z| Vim) = 0. (5) m? a(l + 1)! 
我 们 将 Wig (0) BFE RIA TB AT: Diy = (We ape hee a $ D) 
m 2 : (15) 
Vim (r) = 之 C1 PS, (CD， = re (Vp + vp)? Nnw: 
3 ; (6) AN N 
F 22 (2y) t: |’ i 3 结果 与 讨论 
Prim) = acre e tn "Y nT). 
vz + 1! 31 SRF 


这 里 将 高 斯 参数 v, = Ur 选取 为 一 系列 成 等 比 关系 的 
数 : 


(n=1~n™*), (7) 
最 后 ， 利 用 Rayleigh-Ritz 变 分 法 求解 如 下 方程 : 


max 
n 


> [T, ant + Vint 了 EN, n| cwi =0. (8) 
n'=1 
, Dirac 方程 利用 矩阵 元 可 以 表示 为 


M max 


2 [Tw + Vn — 


A!=] 
+S nt + = EN wn | Cwi = 0, 
(n = l — Ana ) : 


其 中 开 ， 是 动能 项 的 矩阵 元 ， 人 六， 代表 库伦 势 的 矩阵 元 ， 
Rant 是 动能 相对 论 修 正 项 的 矩阵 元 ，S,， 是 自 旋 轨道 耦 
合 项 的 矩阵 元 ，D,,, 代表 Darwin 项 的 矩阵 元 。 具体 矩阵 
元 表示 为 : 


同样 


t 


(9) 


Vn Vy’ 


2 1+4 
G G n 
N, n 一 (Prim | Porm) = ( v FV ) ? (10) 


= G 
Tan = ( nim fe) 


-2 y 


11 
T h2 (21 + 3)v,, Vy Wy 
mm wi OO? 
Vani = (nl lp a) 
1 24+1711 (12) 
SSS ee iN 19 
a Vn E Ma Non 
T Prim Sie ban) 
y , (13) 
= 5021 +3)Q1+5 Dr O) N, 


首先 ， 通 过 计算 氧 原子 能 级 ， 并 将 其 计算 结果 同 精 
确 解 进行 对 比 ， 来 验证 高 斯 展开 法 的 精确 性 〈 见 表 1) 。 
对 于 基态 ， 高 斯 展开 法 的 计算 结果 与 精确 解 在 小 数 点 后 
第 12 位 才 开 始 有 差别 ; 而 对 于 2P 态 ， 其 在 小 数 点 后 15 
位 才 有 差别 , 误差 只 有 107 ~ x107" eV. 接 下 来 , 通过 
求解 非 相对 论 极限 下 Dirac 方程 (2) 来 获得 氧 原子 29 
态 和 2P 态 结合 能 ， 并 将 其 结果 与 一 阶 微 扰 的 结果 进 
行 对 比 〈 见 表 2) ， 结 果 表 明 ， 对 于 2S5 55 2P X, 
阶 微 扰 计算 出 的 结果 在 小 数 


完全 对 角 化 给 出 的 结果 与 
点 后 第 9 位 产生 差别 。 
表 2 展示 了 方程 (2) 中 Dirac 方程 修正 项 的 期 望 值 。 
动能 修正 项 R) 对 于 氢 原 子 能 级 的 影响 只 由 波 函 数 贡 
献 ， 所 以 有 ERl2Sin)> Er(2Pip)，Darwin 项 (D) 只 对 角 
动量 为 零 的 s 波 有 影响 ， 使 得 2.S， 能 级 上 移 。 同 时 自 旋 
轨道 耦合 项 (S) 只 对 非 s 态 有 影响 且 与 总 角 动 量 的 值 相 
关 ， 会 使 得 2P, 态 能 级 下 移 。 
本 小 节 的 计算 表明 高 斯 展开 法 可 以 用 来 计算 研究 非 
相对 论 与 相对 论 效 应 下 的 氢 原 子 能 级 ， 能 够 可 靠 的 区 分 
出 各 种 效应 下 的 氢 原 子 能 级 的 微小 变化 。 
表 1 ARF 1S, 28 与 2P 能 级 的 精确 解 与 高 斯 展开 法 的 计算 
结果 。 (原子 单位 制 ) 


精确 值 GEM 计算 结果 
1s -0.5000000000000000 -0.499999999992102 
2S -0.1250000000000000 -0.124999999999033 
2P -0.1250000000000000 -0.125000000000005 


R2 i e g e 2S BS 
Pi 态 的 计算 结果 。 (原子 单位 制 ) 


2S1 2Pip 


非 微 扰 -0.1250020841 -0.1250020814 
一 阶 微 扰 -0.1250020812 -0.1250020812 
R 0.0000054097 0.0000009712 

D 0.0000033285 0 
S 0 -0.0000011100 
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3.2 织 氢 原子 的 兰 姆 位 移 


本 节 对 缘 氨 原子 进行 了 一 系列 的 计算 ,前 文 提 到 , 相 
比 氨 原 子 ， 由 于 缘 子 的 质量 远 远大 于 电子 ， 缪 氨 原 子 系 
统 的 玻 尔 半径 大 约 是 氧 原子 的 11200。 表 3 中 给 出 比 毛 原 
子 玻 尔 半 径 的 计算 结果 为 285.7 fm, 相 比 氧 原子 的 32918 
fm 小 了 约 200 倍 ; 表 3 的 三 四 栏 为 缘 毛 原子 与 氨 原 子 的 
基态 与 第 一 激发 态 的 能 级 ， 可 以 看 出 缪 所 原子 的 束缚 能 
普遍 比 氨 原子 同样 能 级 的 束缚 能 大 200 倍 左右 。 

表 3 只 考虑 库伦 势 的 情况 下 缪 氢 原 子 与 氨 原 子 玻 尔 半径 ， 基 
态 能 量 与 第 一 激发 态 能 量 的 对 比 。 


玻 尔 半 径 基态 [eV] 第 一 激发 态 
[fm] [eV] 
ZUR 52918.74 -13.61 -3.39 
ZAAT 285.71 -2528.85 -632.21 


接 下 来 , 通过 求解 方程 (3) RREGBARTH 2S 
态 和 2P， 态 的 能 级 。 如 果 只 考虑 缪 子 和 质子 间 库 伦 势 ， 
即使 用 方程 (1) 来 求解 缘 毛 原子 ， 所 得 到 的 2S 态 和 
2Pi, 态 的 能 级 在 10-3 meV 的 量 级 也 是 完全 一 样 的 ( 见 
表 4 第 二 列 )， 而 在 考虑 了 相对 论 修正 后 (包含 质子 半径 
以 及 质子 纪 子 约 化 质量 修正 )， 求 解 Breit 哈密 顿 量 得 出 
了 新 的 缘 氢 原子 的 2S SA 2P 态 的 能 级 〈 见 表 4 第 
三 列 ) 。 发 现 2S 态 和 2Pij, 态 的 能 级 仍然 接近 但 是 有 
微小 的 差别 ， 且 分 别 与 只 考虑 库伦 势 时 的 能 级 产生 了 约 
0.023 meV 的 偏 移 。 

表 4 BIRT 2S 和 2P， 态 的 能 级 。 单 位 为 电子 伏特 (eV)。 


库伦 势 Breit 
2S17 -632.2119203 -632.2343889 
oD, -632.2119204 -632.2343674 


尽管 在 考虑 了 相对 论 修正 后 ，2.51,, 态 和 2P 态 的 
能 级 产生 了 大 约 0.02 个 meV 的 差别 ， 但 是 与 实验 值 


的 影响 会 接着 影响 Uehling 势能 的 计算 结果 。 同时 , 与 只 
考虑 库伦 势 时 的 一 阶 微 扰 计算 结果 类 似 , 考虑 Breit 哈密 
顿 量 时 的 一 阶 微 扰 计 算 结 果 同 样 也 比 非 微 扰 的 计算 结果 
低 0.12 meV 左右 。 综 上 ， 使 用 高 斯 展开 法 的 非 微 扰 计算 
以 及 考虑 相对 论 修正 后 兰 姆 位 移 的 计算 结果 ，205.2970 
meV, 更 加 接近 实验 结果 。 图 1 中 直观 的 展示 了 缪 氨 原 子 
在 只 考虑 库伦 势 ， 考 虑 相对 论 效应 以 及 加 上 Uehling 势 
能 之 后 的 基态 与 第 一 激发 态 的 能 级 。 


最 后 ， 我 们 研究 相对 论 效应 对 质子 半径 的 影响 。 在 
引言 中 提 及 了 可 以 通过 实验 测量 缪 氢 原 子 的 兰 姆 位 移 的 
能 级 臂 裂 来 间接 的 获得 质子 的 电荷 半径 。 由 于 理论 计算 
值 同 实验 值 之 间 存 在 5.00 的 差异 ， 为 了 弥补 这 一 差异 ， 
FERN BS AUR F 28 AS 5 2P 态 的 能 级 差 值 
(4(2S1 2P) RAR TE, 28, pF 
2P 态 的 能 级 差 99% AAG ee HB SUR FAY 
贡献 ， 剩 下 部 分 则 由 256 5 2P, 态 的 超 精 细 结 构 劈 
N, 2P 5 2Py 态 的 精细 结构 劈 裂 贡献 ， 以 及 质子 
电荷 半径 相关 项 [15-16,29-30] : 


TN 


= 


2 3 
rp rp 

E, = | —5.2262— + 0.0347 — | meV (16) 
fm? fm? 


所 贡献 。 其 中 29, 2 5 2P 态 的 超 精 细 结 构 劈 裂 以 及 
2Pip 态 与 2P， 态 的 精细 结构 劈 裂 都 是 由 Breit-Pauli 哈 


密 顿 量 所 贡献 。 


为 了 研究 相对 论 效应 下 兰 姆 位 移 的 变化 对 质子 电荷 
半径 的 影响 , 计算 了 两 种 情况 下 的 AS -2P 情况 
一 ， 将 Breit 哈密 顿 量 给 出 的 29, pA 5 2P 态 的 超 精 
细 结 构 辟 裂 与 2P,, 态 与 2P, 态 的 精细 结构 辟 裂 的 单独 
提取 出 ， 再 加 上 仅仅 考虑 库伦 势 时 的 兰 姆 位 移 ; 情况 二 ， 
直接 计算 Breit 哈密 顿 量 加 Uehling 势能 下 的 ASi — 


206.2949 meV! 仍然 不 在 同一 个 量 级 ， 且 理论 上 相对 
论 效应 并 不 会 导致 兰 姆 位 移 。 因 此 在 缪 氢 原 子 中 ， 为 了 
计算 251, 态 与 2Pj, 态 的 能 级 臂 裂 , 我 们 使 用 Uehling 势 
能 的 Yukawa 近似 的 形式 。 首先 , 仅 使 用 库伦 势 和 Uehling 
势能 来 计算 兰 姆 位 移 , 其 结果 在 表 5 第 二 行 , 计算 出 的 兰 
姆 位 移 为 205.2132 meV， 并 与 一 阶 微 扰 的 结果 进行 对 比 
( 表 5 第 三 行 ), 发 现 一 阶 微 扰 计 算出 的 兰 姆 位 移 要 比 全 对 
角 化 计算 的 结果 约 小 0.12 meV。 最 后 同时 考虑 Breit 哈 
密 顿 量 与 Uehling ARE, 全 对 角 化 后 的 计算 结果 在 表 5 第 
四 行 ， 此 时 ， 与 仅仅 考虑 库伦 势 时 相 比 ， 兰 姆 位 移 的 计 
算 值 产生 了 约 0.09 meV 的 偏 移 。 此 偏 移 产 生 的 原因 如 
下 : 尽管 Breit 哈密 顿 量 中 除了 库伦 势 以 外 的 项 对 2.5), 
AM 2P 态 带 来 的 影响 相 比 库伦 势 项 是 很 小 的 ， 但 是 
很 明显 两 种 情况 下 波 函 数 是 不 一 致 的 ， 波 函数 变化 带 来 


2Psp) 值 。 两 种 情况 的 比较 在 表 6 的 第 二 列 展示 。 可 以 
看 出 这 两 种 情况 的 ASi - 2P) 的 差 值 与 表 5 中 第 二 
行 与 第 四 行 的 兰 姆 位 移 的 差 值 几乎 一 样 。 然 后 ， 将 两 种 
情况 下 计算 出 的 4(2S1 - 2B) 值 分 别 与 实验 测量 出 的 
A(2S a —2 Py) {A (206.2949 meV) PO 对比, 不同 的 部 分 
由 式 16 中 的 质子 电荷 半径 相关 项 贡献 ， 即 可 以 计算 出 两 
种 情况 下 质子 的 电荷 半径 ， 具 体 结 果 在 表 6 的 最 后 一 列 ， 
可 以 看 出 相对 论 效应 下 兰 姆 位 移 的 变化 会 导致 质子 的 电 
荷 半 径 产生 约 0.01 fm 的 差别 。 考 虑 到 目前 一 系列 缘 氨 
原子 能 级 实验 给 出 的 质子 电荷 半径 在 0.84 fm 左右 ， 而 
电子 质子 散射 实验 给 出 的 质子 电荷 半径 则 在 0.87 fm 左 
右 。 因 此 ， 相 对 论 效应 与 对 兰 姆 位 移 的 影响 是 需要 考虑 
的 ， 其 有 可 能 为 解决 两 类 实验 中 质子 电荷 半径 差异 的 问 
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表 5 Uehling 势能 对 缪 氢 原 子 2917 态 与 2Pip 态 能 级 的 修正 ， DR =o it Be TE, STAs 电子 伏特 (meV) o 


42512 A2 Pin A2 Pip = A2Si 
库伦 势 -219.8001 -14.5869 205.2132 
库伦 势 (一 阶 微 扰 ) -219.6793 -14.5842 205.0952 
Breit -219.8878 -14.5908 205.2970 
Breit( 一 阶 微 扰 ) -219.7670 -14.5880 205.1790 
实验 结果 206.2949 
Coulomb Breit Breit + Uehling 
632.18 


E(eV) 


-632.46 


Ds 
a 


Ro BART 25, 55 2P,, SHARE (单位 : 毫 电子 伏 
特 ) 以 及 对 应 的 质子 电荷 半径 x,( 单 位 :fm)。 


A(2.S 1) — 2P;) ry 
情况 一 209.9887 0.8431 
情况 二 210.0778 0.8532 


4 结论 


本 文 首次 利用 高 斯 展开 法 计算 了 相对 论 框架 下 的 缪 
氨 原 子 部 分 能 级 以 及 缪 所 原子 的 兰 姆 位 移 ， 对 于 氧 原子 
和 缪 所 原子 的 相对 论 修正 ， 分 别 使 用 了 非 相 对 论 极 阶 
的 狄 拉克 方程 与 Breit-Pauli FE. MAART N= 
移 用 描述 了 真空 极 化 的 Uehling 势能 作为 其 等 效 势 。 
首先 ， 对 于 氧 原子 系统 ， 本 文 分 别 求解 了 薛 定 廖 
程 与 非 相 对 论 极 限 下 的 狄 拉克 方程 ， 并 将 薛 定 订 方 程 的 


BAUR FMA ERE 


方程 对 28. A 


态 与 2Pi， 态 造 成 


同时 , 与 只 考 


EA 


了 大 约 0.02 meV 的 差异 ; 


十 伦 势 的 情况 对 比 之 后 发 ZEIL, 其 对 2.917 


SI 2P 态 本 身 都 产生 了 约 0.023 eV 的 差别 。 
TERES 77 FEN Uehling 势能 与 Breit- 
Pauli 方程 加 Uehling 势能 下 的 兰 姆 位 移 ， 结 果 表 明 后 


最 后 求解 


果 低 0.12 meV， 即 更 偏离 实验 


者 计算 出 的 兰 姆 位 移 比 前 者 高 了 约 0.09 meV, Hh BRE 
实验 结果 。 另 一 方面 ， 如 果 把 Uehling 势能 看 成 微 扰 项 


计算 一 阶 微 扰 下 兰 姆 位 移 ， 其 结果 比 非 微 扰 的 计算 结 


结果 。 在 此 基础 上 ， 利 用 


公式 16 计 算 了 两 种 情况 下 的 质子 电荷 半径 ， 发 现 相 对 论 


修正 导致 的 这 一 0.09 meV 的 1 


l 移 会 对 质子 电荷 半径 带 


来 约 0.01 fm 的 变化 ， 此 变化 在 目前 一 系列 质子 电荷 半 
径 的 测量 实验 中 是 不 可 忽视 的 。 因 此 ， 在 兰 姆 位 移 的 这 


= 


计算 结果 与 精确 解 进行 了 对 比 、 将 非 相 对 论 极 限 下 狄 拉 


殉 方 程 的 计算 结果 与 一 阶 微 扰 的 结果 进行 对 比 ， 结 果 表 


明 高 斯 展开 法 的 精确 度 在 10-?eV 量 级 。 
接 下 来 ， 分 别 通过 求解 薛 定 谓 方 程 与 Breit-Pauli 方 
FERIA AUR SAE, SHARES IEEE, Breit 


样 对 精确 度 要 求 很 高 的 计算 中 ， 不 宜 使 用 一 阶 微 扰 来 进 
行 计 算 。 同 时 在 纱 氧 原子 中 ， 考 虑 相对 论 修正 是 必要 的 ， 


相对 论 修 1 


MB SURF A 


对 波 函 数 的 影响 进而 会 对 兰 姆 位 移 产生 影响 ， 且 考虑 


正 后 的 计算 结果 更 接近 实验 值 ， 此 结论 对 后 续 
目 关 计算 有 重要 的 参考 价值 。 此 外 ， 由 于 


本 文 只 考虑 J 单 圈 真 空 极 化 对 了 


F 兰 姆 位 移 的 贡献 ， 但 是 
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“5. 


LAHAT, Bo 


9 能 修正 等 (346 都 会 对 兰 姆 


位 移 产 生 大 于 0.09 meV 的 影响 ， 如 果 利 
程 加 上 这 些 修正 项 进行 计算 ， 


用 Breit-Pauli 方 
也 许 会 抵消 或 者 增 大 这 一 


3 


0.09 meV 的 差异 ， 进 而 对 质子 电荷 半径 产生 影响 。 后 续 


会 基于 此 进行 } 


bs 
a 


多 研究 计算 。 高 精度 测量 是 目前 粒 


子 物 理学 重要 的 实验 前 沿 之 一 ， 本 文 对 缪 氢 原 子 兰 姆 位 


移 的 计算 将 为 缪 子 物理 的 高 精度 测量 实验 提供 重要 的 班 
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High precision computation of the Lamb shift of the muonic hydrogen atom 


within the relativistic framework 
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Abstract: The measurement of the muonic hydrogen atom Lamb shift is one of the most fundamental ways to obtain 
the proton charge radius experimentally. The theoretical calculation of the muonic hydrogen atom Lamb shift is nec- 
essary to extract the proton charge radius from the Lamb shift. In this work, we calculate the muonic hydrogen atom 
Lamb shift within the relativistic framework. The relativistic effect is described with the Breit-Pauli hamiltonian. The 
vacuum polarization effect is realized with the effective Uehling potential. With a high precision calculational method, 
Gaussian expansion method, we calculate the 291， and 2P levels of the muonic hydrogen atom. We find that the 
muonic hydrogen atom Lamb shift with considering the relativistic effect is about 0.09 meV larger than the one without 
considering the relativistic effect and the former one is closer to the experimental value. 


Key words: muonic hydrogen; lamb shift; Gaussian expansion method 
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